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LARGE-EDDY  SIMULATION  OF  MACH  3.0 
FLOW  PAST  A 24-DEGREE  COMPRESSION  RAMP 


DONALD  P.  RIZZETTA  AND  MIGUEL  R.  VISBAL 
Air  Force  Research  Laboratory 
Wright-Patterson  Air  Force  Base,  Ohio  45433-7521 


Abstract. 

Large-eddy  simulation  of  the  Mach  3.0  flow  past  a 24  deg  compression 
ramp  is  performed  by  a high-order  numerical  method.  Spatial  derivatives 
are  represented  by  a sixth-order  compact  stencil  that  is  used  in  conjunction 
with  a tenth-order  non-dispersive  filter.  The  scheme  employs  a time-implicit 
approximately-factored  finite-difference  algorithm,  and  applies  Newton-like 
subiterations  to  achieve  second-order  temporal  and  sixth-order  spatial  ac- 
curacy. In  the  region  of  the  shock  wave,  compact  differencing  of  convective 
fluxes  is  replaced  locally  by  a third-order  Roe  upwind-biased  evaluation. 
The  Smagorinsky  dynamic  subgrid-scale  model  is  incorporated  in  the  sim- 
ulation to  account  for  the  spatially  under-resolved  stresses  and  heat  flux. 
Comparisons  are  made  with  experimental  data  in  terms  time-mean  sur- 
face pressure  and  skin  friction  distributions,  and  with  instantaneous  surface 
pressure  measurements. 


1.  Introduction 

Large-eddy  simulation(LES)  of  supersonic  flows  is  useful  for  studying  com- 
pressibility effects,  which  can  appreciably  alter  fluid  physics.  Such  studies 
increase  the  understanding  of  turbulence  mechanisms  and  can  lead  to  the 
development,  improvement,  and  testing  of  lower-order  closure  models.  De- 
spite remaining  computationally  intensive,  LES  also  may  be  beneficial  in 
the  design  and  analysis  of  high-speed  flight  vehicles  and  associated  propul- 
sion systems  where  less  sophisticated  approaches  fail. 

Due  to  their  geometric  simplicity,  supersonic  compression-ramp  flow- 
fields  have  been  studied  extensively,  both  experimentally  and  computa- 
tionally. Characteristics  of  the  unsteady  shock-wave  motion  of  such  flows 
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have  been  observed,  measured,  and  analyzed  by  Andreopoulos  and  Muckfl], 
Smits  and  Muck[2],  Dolling  and  Murphy[3],  and  Erengil  and  Dolling[4, 
5],  among  others.  Numerical  investigations  have  typically  considered  the 
Reynolds-averaged  Navier-Stokes  equations  used  in  conjunction  with  mean 
turbulence  models.  These  efforts  have  met  with  limited  success  in  the  pre- 
diction of  quantities  such  as  heat  transfer  and  skin  friction,  particularly  in 
situations  with  large  reversed-flow  regions[6].  It  is  believed  that  this  dif- 
ficulty may  be  due  in  part  to  the  disparity  between  the  time-mean  and 
instantaneous  shock-system  structure.  In  addition,  the  models  and  resul- 
tant computations  often  fail  to  account  for  compressibility  effects  or  the 
three-dimensionality  of  the  flowfield.  In  an  effort  to  overcome  these  defi- 
ciencies, direct  numerical  and  large-eddy  simulations  have  been  carried  out 
by  Hunt  and  Nixon[7],  Urbin  et  al.[8,  9],  Adams[10,  11],  and  Rizzetta  et 
al.[12,  13]  The  present  effort  provides  a large-eddy  simulation  for  the  flow 
past  at  24  deg  compression  ramp  at  a freestream  Mach  number  of  3.0  and 
Ree  = 1696. 


2.  Governing  Equations 

The  governing  equations  are  the  unsteady  three-dimensional  compressible 
Favre  filtered  Navier-Stokes  equations,  written  in  nondimensional  variables 
utilizing  a generalized  coordinate  system,  and  expressed  notationally  in  the 
following  conservative  form 


§ + |<F  - 5* ’•>  + - ib'* + - jb*>  = M 


All  dependent  variables  may  be  decomposed  into  a filtered  or  large-scale 
portion,  and  a subgrid-scale  component  e.g.,  u = u + usg.  It  is  then  con- 
venient for  compressible  flows  to  recast  the  large-scale  component  in  terms 
of  a Favre-averaged  variable  so  that 


pu 

u = — . 

P 

The  subgrid-scale  stress  and  heat  flux  are  provided  by 

Tij  = -Rep(uiUj  - UiUj),  Qi  = Rep(uiT  - UiT). 


(2) 

(3) 


Following  Germano  et  al.[14],  the  compressible  version  of  the  model  is 
given  in  trace-free  form  as 


1/2 


pt  — ReCA2p§M  where  Sm  — 


(4) 
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is  the  magnitude  of  the  rate-of-strain  tensor,  and 

c..  = 1 ( , djk  duj\ 

2 \dxjd£k  dxidtk) 


(5) 


The  eddy-viscosity  length  scale  is  taken  as 


which  corresponds  to  the  width  of  the  grid  filter  in  physical  space,  C is  the 
eddy-viscosity  model  coefficient,  and 


For  compressible  applications,  the  isotropic  part  of  the  stress  tensor  is  ob- 
tained according  to  Yoshizawa[15]  from 


rkk  = 2CIA2pS2M.  (8) 

To  complete  closure  of  the  model,  the  subgrid-scale  heat  flux  vector  is 
specified  in  terms  of  the  turbulent  Prandtl  number  as 


\Pr,J  dx,d(j 

The  model  coefficients  C,  Ci , and  the  turbulent  Prandtl  number  Prt  are 
computed  as  a function  of  time  and  space  from  the  energy  content  of  the 
resolved  large-scale  structures. 


3.  Numerical  Method 

Time-accurate  solutions  to  Eq.  1 were  obtained  numerically  by  the  implicit 
approximately- factored  finite-difference  algorithm  of  Beam  and  Warming[16] 
employing  Newton-like  subiterations,  which  may  be  represented  notation- 
ally  as  follows 
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-(T)  [(4)(3°p-4Q"+«”_1) 
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In  this  expression,  which  was  employed  to  advance  the  solution  in  time, 
Qp+1  is  the  p+l  approximation  to  Q at  the  n+1  time  level  Qn+1,  and  A Q = 
Qp+l  — Qp.  The  implicit  segment  of  the  algorithm  incorporated  second- 
order-accurate  centered  differencing  for  all  spatial  derivatives^^,  <^2,  <^2), 
and  utilized  the  diagonalized  form  of  the  factorized  equations  to  enhanced 
efficiency.  Temporal  and  spatial  accuracy  were  maintained  by  utilizing  subit- 
erations within  a time  step. 

A sixth-order  tridiagonal  subset  of  Lele’s[17]  compact  difference  scheme 
was  used  to  evaluate  the  right-hand  side  of  Eq.  10(^6,  <^6,  <^6),  and  is 
illustrated  here  in  one  spatial  dimension  as 


'dF\ 

- +a 


= a 


<%lAsa 


dF\ 

dUm 


(11) 


with  a = 1/3,  a = 14/9,  b — 1/9.  (12) 

The  scheme  is  used  in  conjunction  with  a lOth-order  non-dispersive  compact 
spatial  filter  in  order  to  maintain  both  stability  and  accuracy,  particularly 
on  stretched  curvilinear  meshes.  The  filter  is  applied  to  the  solution  vector 
sequentially  in  each  of  the  three  computational  directions  following  each 
subiteration,  and  is  implemented  as 


&fQi- 1 + Qi  + OlfQi+ 1 — ^2  A(Qi+n  + Qi-n)  (13) 

71  = 0 L 

where  Q is  the  filtered  value  of  Q.  Equation  13  represents  a one-parameter 
family  of  filters,  where  numerical  values  for  a/  and  the  an’s  may  be  found 
in  Ref. [18].  Repeated  application  of  the  spatial  filter  can  result  in  shock 
waves  that  are  excessively  diffuse.  This  deficiency  is  overcome  by  replacing 
the  compact-differencing  of  convective  derivatives  and  use  of  filtering,  by 
Roe’s  third-order  upwind-biased  scheme,  locally  in  inviscid  regions  of  shock 
waves. 
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4.  Results 

A computational  domain  size  was  taken  as 

Lx  = 31.2(50,  Ly  = 4.7<50,  Lz  = 2.950,  (14) 

where  Lx,Ly,Lz,  are  streamwise,  vertical,  and  spanwise  extents  respec- 
tively, and  Sq  is  the  height  of  the  incoming  boundary  layer.  The  vertical 
extent  Ly  corresponds  to  the  domain  height  at  the  inflow  location.  This 
region  was  discretized  with  a nonuniform  computational  mesh  consisting  of 
(421  x 151  x 81)  points  in  (i,j,k).  At  the  inflow  location,  the  grid  had  the 
following  minimum  spacings  in  wall  units: 

Ax+  = 16.8,  A y+  = 1.4,  Az+  = 8.3.  (15) 

Based  upon  the  mean  incoming  profile,  79  of  the  151  vertical  grid  points 
were  within  the  boundary  layer. 

The  solution  presented  here  for  the  Mach  3.0  flow  past  a 24  deg  com- 
pression ramp  with  Reg  = 1696  represents  part  of  a comprehensive  investi- 
gation which  was  performed  for  supersonic  compression-ramp  flows.  These 
are  described,  along  with  more  complete  details  of  the  present  simulation, 
in  Ref.  [13].  Typical  instantaneous  results  at  the  midspan  location  are  pre- 
sented in  Fig.  1.  Mach  number  contours  appear  in  the  left-hand  portion  of 
the  figure,  while  the  grid  structure  is  indicated  to  the  right.  Segments  of 
the  grid  which  are  blanked  out  correspond  to  mesh  points  where  convective 
derivatives  were  obtained  via  the  Roe  upwind-biased  scheme.  In  the  figure, 
only  every  other  f-grid  and  every  third  j-grid  line  are  displayed.  It  is  evi- 
dent that  the  upwind-biased  evaluation  is  confined  locally  to  a small  region 
surrounding  the  shock  wave,  so  that  accuracy  of  the  high-order  method  is 
not  compromised  in  other  regions  of  the  flowfield. 

Although  no  experimental  data  exists  for  the  exact  flow  conditions  of  the 
simulation,  it  is  useful  to  compare  to  measurements  obtained  at  Reynolds 
numbers  which  were  several  orders  of  magnitude  larger  than  that  of  the 
computation.  Spanwise  averaged  time-mean  surface  pressure  and  skin  fric- 
tion coefficient  distributions  are  shown  in  Fig.  2.  The  pressure  coefficient 
has  been  normalized  by  the  inviscid  rise  to  account  for  variations  in  the 
freestream  Mach  number  between  the  calculation  and  experiments.  Be- 
cause of  the  higher  Reynolds  numbers  of  the  measurements,  Cf  has  been 
normalized  by  its  value  just  upstream  of  the  interaction  region. 

Comparisons  of  the  computed  surface  pressure  standard  deviation(.s) 
and  skewness(5fc)  distributions  with  the  data  of  Dolling  and  Murphy[3] 
appear  in  Fig.  3.  The  value  of  s upstream  of  the  interactional ) has  been 
removed  from  the  standard  deviation  in  order  to  account  for  differences  in 
the  incoming  states  between  the  simulation  and  the  experiment.  Apart  from 
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disparities  near  separation(X/^o  ~ —2.0),  these  comparisons  are  favorable. 
The  disparities  are  caused  by  differences  in  the  shock  wave  motion  between 
the  respective  results,  which  will  be  illustrated  subsequently. 

The  wall  pressure  intermittency(T)  and  instantaneous  pressure  time 
history  are  displayed  in  Fig.  4.  It  is  observed  that  the  experimental  dis- 
tributions of  F have  a steeper  rise  which  occurs  further  downstream  than 
the  numerical  result.  This  is  due  to  the  more  extensive  interaction  region 
present  in  the  low  Reynolds  number  large-eddy  simulation.  Although  the 
time  mean  level  of  the  computed  surface  pressure  near  the  separation  point 
is  approximately  the  same  as  that  of  the  experiment,  a very  different  fluc- 
tuating component  is  indicated.  The  high  frequency  oscillations  exhibited 
in  the  numerical  simulation  are  similar  to  those  of  the  incoming  boundary 
layer,  while  the  low  frequency  modes  of  the  experiment  correspond  to  more 
extensive  motion  of  the  shock  wave. 


5.  Summary 

Computed  surface  pressure  and  skin  friction  distributions  compared  rea- 
sonably well  with  experimental  data  collected  at  higher  Reynolds  numbers. 
Comparisons  were  also  made  with  statistical  quantities  extracted  from  in- 
stantaneous unsteady  surface  pressure  measurements. 
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Figure  1.  Typical  instantaneous  Mach  number  contours  and  shock-capturing  stencils 
at  the  midspan  location 
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Figure  2.  Spanwise  averaged  time-mean  surface  pressure  and  skin  friction  coefficient 
distributions 


Figure  3.  Spanwise  averaged  surface  pressure  standard  deviation  and  skewness  distri- 
butions 


Figure  4 ■ Spanwise  averaged  surface  pressure  intermittency  distributions  and  midspan 
surface  pressure  time  history  at  X/So  = —2.1 


